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Abstract 

I predict the existence of internal spatial currents in a single macroscopic 
quantum system, namely in trapped dilute-gas at sufficiently low tempera- 
tures, when a Bose-Einstein condensation occurs. The spatial profiles of the 
wavefunctions of low-lying states in such a system are different due to the 
inhomogeneity, caused by an asymmetry of external trapping potential. This 
is the reason for appearing of Josephson-like oscillations between atomic sub- 
systems in different states including the ground state as well. Using a simple 
model for the wavefunctions of three low-lying states we demonstrate how 
essential this effect can be. The possible applications of the predicted effect 
are briefly discussed. Particularly, this effect opens the possibility to identify 

experimentally the low lying excited states of a system. 
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The recent observations of Bose-Einstein condensation (BEC) in trapped atomic gases 
|jlj] have renewed interest in bosonic many-particle systems at low temperatures. The strong 
inhomogeneity of trapped gases in real experiments makes a system to be even more inter- 
esting. 

Of course, the most spectacular consequences of Bose-Einstein condensation is superflu- 
idity and high phase coherence 0. The latter can lead to quite a number of interference 
phenomena, which were not being observed earlier in macroscopic systems in such a direct 
way. One of such interference effects was observed by the MIT group || . They used a laser 
beam to separate an atomic cloud into two parts. After switching off the confining poten- 
tial and the laser, the authors of || observed clean interference patterns in the overlapping 
region. 

The authors of showed a possibility of another interesting manifestation of the phase 
coherence, namely they predicted the existence of Josephson-type effects in a double-well 
external potential. In this case a difference of chemical potentials ,say \L\ and \x-i-, in two 
condensates leads to appearance of the Josephson current / oc sin[(/ii — H2)/hi\. Analytic 
calculations of the spatial current distribution are rather complicated because of necessity 
to account the boundary conditions of condensates order parameters 0, ||. Nevertheless, 
the Josephson oscillations are already observed in experiment || and are extensively studied 
by various groups |J both theoretically and experimentally. 

Here we propose a simple model theory for internal spatial oscillations, which occur in a 
single condensate. The origin of this effect is a difference of spatial profiles of wavefunctions 
of the ground state and neighboring excited states. So, the well-known quantum interference 
between two quantum states can lead to the observable spatial process in the systems under 
consideration. 

It was realized long ago, that the ground state properties of trapped Bose-gas can be 
well described by the Gross-Pitaevskii equation J7J 

ih ^m^ = {~H A + Utrap{r) + 9Mri t)|2 } $(r ' t)} (1) 
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where $(r, t) = (ip(r,t)) is the condensate wavefunction, Ut Ta . P (r) is the trapping potential, 

4:7ih 2 a 

9 = 

m 

is the zero-range interaction constant, a is the s-wave scattering length. To obtain the ground 
state properties, one can write the condensate wavefunction as $(r, t) = ip(r) exp(— ifxt/h). 
In fact, this mean we assumed a condensate to be in the stationary state. In this case 
equation ([I]) takes the stationary form 

Mr) = l~A + C/ trap (r) + ^(r)| 2 j <p( r ). (2) 

This equation explains the sense of the chemical potential fi as an energy of the stationary 
level. This can be used further for determining of ip(r) and /i. 

In this Letter we want to study the ground state of trapped bosons in an asymmetric 
harmonic potential. For our aims we must allow for (at least) two single-particle states. 
Keeping in mind the illustrational purposes we shall consider the simplest case of stationary 
states, where the single-particle wavefunctions have the oscillator-like form || 

^(r,*) = ^(r)exp(-i^), (3) 

where 

^o(r) = -^exp(-^), (4) 
V2 ( r 2 \ 

b is a variational parameter. These have the same shape as the ground state and the first 
excited state with a zero projection of angular momentum. More generally, we could take 
</9i(r) as a linear combination of the states with different angular momentum projections 
mi = 0, ±1, but this has no influence on the qualitative picture. Furthermore, below we 
will operate in fact with the one dimension, say z. This is particularly acquitted in view of 
spatial anisotropy of trapping potential in majority of real experiments 0. The resulting 
wavefunction of our simple model system can be presented as a linear combination: 



& 01 \v,t) = ^o(r)exp (-i^) + ^(r) exp (-i^tj . (6) 

The current density is then given by the well-known relation [T(J 

l(oi) = th_ {$(oi)( r)t ) V $(w)*(r,t) - $( 01 )*(r,t)V$ (01) (r,t)} . (7) 

Substituting eq. (ffl) into eq. (0), we come to 

I^ = --^o(r)^(r)^sin(^^t), (8) 
m z \ n / 

where e 2 is the unit vector along z-direction. Equation (|j) looks very similar to the cor- 
responding expression in |4[] for the Josephson current between two condensates, but have 
rather different sense. This simple example shows that there is an oscillating current which 
results in oscillations of particle number in each state. For example the temporal variation 
of the number of ground state particles due to a Josephson-type exchange with the first 
excited state is given by the formula 

Typical dependence of Sn^\z,t) on z and t is presented in Fig. [I]. 

Significance of the discussed effect is especially clearly manifests if we account the second 
excited state as well. Taking the simple oscillator-type wavefunction [TIJ for this second 
excited state, we easily get the following formula 

5ny{z,t) = -exp -— cos 1 , (10) 



7T 3 / 2 m& 6 (yL*2 — A*o) \ b 2 J V h 
where \x stands for the chemical potential of the second excited level. Fig. |2| shows the 
dependence of a total variation of the particle number in the ground state 6no(z, t) = 
8no(z, t) + 8riQ(z, t) on time (in units of h/ — /io)), calculated from expressions @ and 
(p!0|). The figure clearly illustrates possible applications of the predicted effect. Namely, if 
we measure the oscillations of the number of particles in the ground state (at any convenient 
coordinate), then we get a possibility to recover the low lying excited stationary states (i.e. 



to find their chemical potentials, spreading, spatial shape etc.). This becomes even more 
tempting in view of a difference in spatial dependences of dn^ and Sn^, which is shown on 
Fig. |3]. Thus, all excited states become well distinguishable by measuring the amplitude and 
period of oscillations. Fig. ^ represents 5uq as a function of both z and t. So, the character 
of oscillations is rather different for each mode. 

In conclusion, using a simple model for describing the ground and low-lying excited states 
of the trapped Bose-gas, we showed the principal possibility of the spatial oscillations within 
a system due to the quantum interference between states with different spatial profiles. These 
oscillations demonstrate rather different behaviour for different energy levels of the system. 
This fact provides a tempting possibility to determine the spectrum of stationary states in 
the system. 

In reality the wavefunctions are of course different from those used above due to an 
interaction. Furthermore, we even do not certainly know whether or not the stationary 
excited states with well defined chemical potentials exist. The measuring of the considered 
internal oscillations just can answer this question. 
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FIGURES 

FIG. 1. The dependence of on z and t (here r = (fii — /j,o)t/fi). We used the constants 
b = 5 ■ 10~ 3 cm, fiQ = 830/i erg, fi\ = 1.3 • 10 3 h erg. Note, the real amplitude of oscillations may 
be, of course, not so large. The parameters of wavefunctions (||) are taken rather approximately 
to draw the general picture. 

FIG. 2. Temporal dependence of 5no = Sn^ + 5n^ at z = b/2 (r = (fi\ — fio)t/h). The 
constants are 6 = 5- 10 -3 cm, fig = 830k erg, [i\ = 1.3 • 10 3 ?i erg, [12 = 3.6 • 10 3 ?i erg. 

FIG. 3. Spatial dependence of 5n^ and Sn^ at r = — /io)t/h = 6. All constants are the 
same as in Figs. [I] and ^. 

FIG. 4. The dependence of 5uq = Sn^ + Sn^ on z and t. All constants are the same as in 
Figs. and §. 
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